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Abstract. Postburst evolution of adiabatic fireballs that 
produce 7-ray bursts is studied. Emphasis has been put 
on the transition from the highly relativistic phase to the 
non-relativistic phase, which, according to our calculation, 
should happen much earlier than previously expected. The 
theoretical light curves decline a little sharper at the non- 
relativistic stage than at early times, but still can fit 
the observations well. However, disagreements are obvi- 
ous when t > 100 d, implying a large initial energy Eq, 
a low interstellar density n, or the possible existence of a 
persistent energy source at the center of the fireball. 
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1. Introduction 

Since their discovery about thirty years ago (Klebesadel 
et al. 1973), 7-ray bursts (GRBs) have made one of the 
biggest mysteries in astrophysics (Fishman & Meegan 
1995). The Italian-Dutch BeppoSAX satellite opened up 
a new era in the field in early 1997. By the end of May 
1998, afterglows have been observed in X-rays from about 
a dozen events, in optical wavelengths in several cases 
(GRB 970228, 970508, 971214, 980326, 980329, 980519, 
and possibly GRB 980425), and even in radio from GRB 
970508, 980329, and GRB 980519. Very recently, GRB 
971214 was reported to have released an enormous energy 
of 3 x 10 53 erg in 7-rays alone (Kulkarni et al. 1998; Wi- 
jers 1998); and GRB 980425 seems to be associated with a 
supernova (Galama et al. 1998b), these facts have aroused 
researcher's more fever and may provide crucial clues to 
our understanding of GRBs. 

The fireball model has become the most popular and 
successful model of GRBs (Meszaros et al. 1994; Fenimore 
et al. 1996, and references therein). After producing the 



main GRB, the fireball will continue to expand as a thin 
shell into the interstellar medium (ISM), generating an 
ultra-relativistic shock. Afterglows at longer wavelengths 
are produced by the shocked ISM (Meszaros & Rees 1997; 
Waxman 1997a,b; Tavani 1997; Sari 1997; Wijers et al. 
1997; Huang et al. 1998; Dai & Lu 1998a,b). It has been 
derived that for adiabatic expansion, R cx i 1 / 4 , 7 oc i -3 / 8 , 
where R is the shock radius measured in the burster's 
static frame, 7 is the Lorentz factor of the shocked ISM 
measured by the observer and t is the observed time. These 
scaling laws are valid only at the ultra-relativistic stage. 

The purpose of this Letter is to study numerically 
the full evolution of adiabatic fireballs from the ultra- 
relativistic phase to the non-relativistic phase. It is found 
that radiation during the mildly relativistic phase (2 < 
7 < 5) and the non-relativistic phase (7 < 2), which was 
obviously neglected in previous studies, is of great impor- 
tance. 

2. Previous studies 

In case of adiabatic expansion, the shocked ISM's Lorentz 
factor evolves based on (Sari 1997; Waxman 1997a, b; Ta- 
vani 1997): 
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where Eq = E51 x 10 erg is the original fireball en- 
ergy, of which about a half is believed to have been re- 
leased as 7-rays during the GRB phase, n = uq 1 cm~ 3 
is the number density of the unshocked ISM, and t s is 
t in units of second. R(t) can be derived from 7 2 i? 3 w 
Eq 1 1 (AimrripC 2 ) , where m p is the proton mass and c the 
velocity of light. Flux density at observing frequency v 
then declines as S v oc < 3 ( 1_ p)/ 4 , where p is the index char- 
acterizing the power-law distribution of the shocked ISM 
electrons, dn' e /d^ c oc 7,7 P . These expressions are valid only 
when 7 > 1. In general, X-ray and optical afterglows 
were observed to follow power-law decays, and such a fire- 
ball/blastwave model agrees with observations quite well. 
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Fig. 1. Evolution of the fireball's Lorentz factor 



However, we notice that afterglows from GRB 970228 
and GRB 970508 have followed simple power-law decays 
for as long as 190 days and 80 days respectively, while 
in Eq.(l), even t = 30 d will lead to 7 <~ 1. We stress 
that the overall evolution of the postburst fireball can not 
be regarded as a simple one-phase process. One should 
be careful in applying those scaling laws at later times. In 
fact, it is clear from Eq.(l) that the expansion will become 
mildly rclativistic (2 < 7 < 5) when t > 10 d, and will 
cease to be relativistic about 30 days later at the latest. 

3. Our Model 

We now propose a refined model to describe the full evo- 
lution of the postburst fireballs. As usual, we suppose that 
the main GRB occurs at, 
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where 77 = 300t7 3O o = E /(M c 2 ), M is the mass of 
the contaminating baryons. In the subsequent expansion, 
jump conditions for the shock can be described as (Bland- 
ford & McKee 1976): 



77 



-n. 



7-1 



77 + 1 , -.x 2 

-(7- l)nm p c , 

7—1 



(7+l)[ 7 (7-l) + l] s 
7(2-7)(7-l) + 2 



(3) 



(4) 



(5) 



where n' and e' are the electron number density and en- 
ergy density of the shocked ISM respectively in the frame 
co-moving with the shell, T is the Lorentz factor of the 
shock, and 7 is the adiabatic index of the ISM, for which 
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Fig. 2. Evolution of the shock's velocity. Full line is our result 
and dash-dotted line is analytic solution for non-relativistic 
shock wave (Eq.(9)) 
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Fig. 3. The fireball's radius vs. time. Full line is our result and 
dash-dotted line is analytic solution (Eq.(8)) 



we have derived an approximation, 7 s=a (47+l)/(37), con- 
sistent with the requirement that 7 ~ 4/3 for an extremely 
relativistic blastwave and 7 i=a 5/3 for the non-relativistic 
Sedov shock. For 7 3> 1, Eqs.(3) — (5) reduce to n 1 = Ajn, 
e' = A"f 2 nm-pC 2 and T = V27, which are just the starting 
point of previous studies. Here we expect that our equa- 
tions are valid for describing relativistic shocks as well as 
non-relativistic blastwaves. 



The kinetic energy of the shocked ISM in the fireball is 
E k = al3 2 T 2 (4/3)irR 3 nm p c 2 (Blandford & McKee 1976), 
where (3 = (1 — 1/r 2 ) 1 / 2 , and a is a coefficient: a — > 0.35 
when (3 — > 1 and a — > 0.73 when (3 — ► 0. We will use an 
approximate equation for a: a = 0.73 — 0.38/3. As usual, 
the expansion is assumed to be adiabatic, during which 
energy is conserved, so we have, 



^TTaP 2 T 2 R 3 nm p c 2 = E k = ^-. 
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Fig. 4. Optical afterglows in R band: a GRB 970228, b GRB 
970508. Full lines are our theoretical results. Observed fluxes 
are from Pedersen et al.(1998), Galama et al. (1998a), Garcia 
et al.(1998) 



In order to study the evolution of 7, we should add the 
following differential equation (Huang et al. 1998), 



dR _ yr 2 - 1 

~dt ~ f 
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Eqs.(2) — (7) present a perfect description of the 
shock. Given initial values of Eq and Mq, R(t) and 7(i) 
can be evaluated numerically. But under the assump- 
tion that 7 S> 1, we can derive a simple analytic so- 
lution, R 4 w R 4 + 8kt, 7 w (cR 3 /k)- 1 / 2 , where k = 
Eq / (Airnm v c) . Additionally if R » i? , that is t » r, 
where r refers to the duration of the main GRB, then R w 
1.06 x lO^^si/no) 1 / 4 * 1 / 4 cm, 7 w 273(£ 51 /n ) 1 / 8 i- 3 / 8 . 
These expressions are consistent with previous studies. 

4. Numerical results 

We have evaluated the propagation of the blastwave 
numerically, taking Eq = 10 51 erg, n = 1 cm -3 , and 
Mq = 2 x 10~ 6 M . Fig. 1 is the evolution of 7, which 
follows the power-law expression of Eq.(l) quite well when 
10 1 s < t < 10 6 s. However, it is clear that 7 ceases to be 
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Fig. 5. X-ray afterglows (0.1 — 10 keV, in unit of erg cm~ 2 
s _1 ). Full line is our theoretical result. For the observed fluxes, 
please see Dai et al.(1998) 



much larger than 1 when t > 10 6 s, so that Eq.(l) is no 
longer applicable. At least, it is problematic to assume 
that the shock wave was still ultra-relativistic, therefore 
one should be cautious in applying those scaling laws such 
as R cx t 1 / 4 , 7 cx t~ 3 / 8 , and S v cx t^ 1 -?)/ 4 . This fact has 
been completely ignored in previous studies. 

Fig. 2 illustrates the evolution of the shock wave's ve- 
locity (V s ). V s differs markedly from c after ~ 10 6 s, and 
the evolution enters the non-relativistic phase, for which 
a simple analytic solution is available (Lang 1980): 



R = 1.15(E k t 2 /p)^ 5 , 



V s = 3R/{10t), 



(8) 



(9) 



where p — nm p . The dash-dotted line in Fig. 2 is plotted 
according to Eq.(9). Our result is in good agreement with 
it. Fig. 3 shows our evolution of radius, also plotted are 
results from Eq.(8). We see from Figs. 3 and 4 that our 
refined model is applicable for both relativistic and non- 
relativistic expansion. 

To compare with observations, we need to calculate 
synchrotron radiation from the shocked ISM. As usual, 
electrons in the shocked ISM are assumed to follow a 
power-law distribution, with an index of p; the magnetic 
field energy density in the co-moving frame is supposed to 
be a fraction £| of the energy density, B' 2 /8tt — £|e'; and 
the electrons are supposed to carry a fraction £ e of the 
energy, ~f Cim i n m c c 2 = £ c ~fm p c 2 . The luminosity distance 
to the GRB source is designated as Dl- We have plotted 
in Figs. 4 and 5 the calculated R band flux densities (<Sr) 
and 0.1 — 10 keV fluxes (Fx) respectively, and compared 
them with observed afterglows. 

In Fig. 4a, we take p — 2.5, £| = 0.01, £ = 1, and 
Z?l = 3 Gpc. Other parameters such as Eq, n, Mq are the 
same as in Fig. 1. Generally the theoretical light curve 
fits GRB 970228 well. However, a problem appears at 
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later times (t > 10 7 s) when the blastwave becomes non- 
relativistic. Our model predicts a sharper decline, consis- 
tent with the conclusion that has been drawn by Wijers et 
al. (1997), but the observed flux decays obviously slower. 
An overall power-law decline lasting for more than 10 7 s 
would require a much larger Eq or a much smaller n to 
ensure that the blastwave is uniformly ultra-relativistic. 
However, due to the limited observational data points, it 
is arbitrary to make any further conclusions. 

For the optical afterglow from GRB 970508, a set of 
physical parameters have recently been carefully derived 
as: E = 3.7 x 10 52 erg, n = 0.035 cm- 3 , £ c = 0.13, 
£| = 0.068, p = 2.2, £> L ~ 4 Gpc (Wijers & Galama 1998; 
Galama et al. 1998c, d). In Fig. 4b, we have taken these 
values, with M = 7.4 x 10~ 5 M Q . Since E is large and 
n is small, the period during which the ultra-relativistic 
approximation is valid is substantially extended. In fact, 
7 keeps to be larger than 1.5 even when t > 1.2 x 10 7 
s, and the theoretical light curve are well presented by a 
single straight line, which fits observations satisfactorily. 
We note that if p = 2.3 were assumed, the observational 
data could be reproduced even better. The afterglow from 
GRB 970508 seemed to decay more slowly after about 
80 days, implying the presence of a constant component 
(Pedersen et al. 1998; Galama et al. 1998a; Garcia et al. 
1998), maybe its host galaxy. The observed optical flux 
peaked about two days later, possibly associated with an 
X-ray outburst (Piro et al. 1998). This feature could not 
be explained by a simple blastwave model. Our theoretical 
light curve peaks several hours later at optical wavelengths 
and several tens of seconds later in X-rays. 

The theoretical X-ray light curve is plotted and com- 
pared with observations in Fig. 5, where parameters are 
evaluated the same as in Fig. 4a. X-ray counterparts usu- 
ally fade away in less than two weeks when the expansion 
is still relativistic, so observed light curves are well repre- 
sented by power-law relations. 

5. Discussion 

The simple fireball model predicts a power-law decay for 
GRB afterglows, in rough agreement with observations. 
However, the fact that the blastwave might enter the non- 
relativistic phase when the afterglows are still detectable 
has been widely ignored. Here we stress that the blast- 
wave generated by a typical fireball will cease to be ultra- 
relativistic in ~ 10 6 s. One should be cautious in applying 
those simple scaling laws such as 7 oc i~ 3 / 8 , R cx i 1 / 4 , and 
S v oc t 3 ' 1 ^)/ 4 (for adiabatic expansion) at such late times 
as t > 10 6 s. 

We have derived refined equations to describe the 
full evolution of postburst fireballs. Our model is consis- 
tent with previous studies in both the ultra-relativistic 
phase and the non-relativistic phase. The predicted after- 
glows decay slightly sharper at the non-relativistic stage. 
The observed overall power-law decay lasting for several 



months is usually considered as a strong proof to the pop- 
ular fireball model. Here we point out that it has in fact 
raised a problem, since it requires a large Eq and/or a low 
n. For example, recent detection of an enormous burst 
GRB 971214 with £ ~ 3 x 10 53 erg might ameliorate 
the ultra-relativistic assumption to some extent (Kulka- 
rni et al. 1998; Wijers 1998). The solution could also be 
that there is a persistent energy source at the center of 
the GRB source (Dai & Lu 1998b). 
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